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Abstract 

The results of perturbative QCD evaluation of the ~ rnj/M'^ contributions to T 
and T z^hadrons for the quark masses ruf <^ Mz are presented. The recent results 
due to the combination of renormalization group constraints and the results of several 
other calculations are independently confirmed by the direct computation. Some existing 
confusion in the literature is clarified. In addition, the calculated 0{a'^) correction to the 
correlation function in the axial channel is a necessary ingredient for the yet uncalculated 
axial part of the 0(af ) mass correction to the Z decay rates. The results can be applied 
to the r hadronic width. 



Recent analyses show that the LEP result for the branching fraction T ^^^^i^/T z^hadrons = 
0.2208 ± 0.0024 ]I[ differs (is larger than) the Standard Model prediction § by 2-2.5 a with 
the top mass at around 174 GeV 0. Although the search for implications of this fact beyond 
the Standard Model is already began (see, e.g., 0), a further analysis within the Standard 
Model is still an important issue. 

Briefly, a current state of the perturbative QCD evaluation of the T z^hadrons and related 
quantities is as follows. The QCD contributions are evaluated to 0{al) in the limits ruf = 
(/ = u,d,s,c,b) and ^ oo p. The leading correction to the above results due to the 
0{c(l) triangle anomaly type diagrams with the virtual top quark has also been calculated in 
the limit ^ oo 0. The calculations to 0(0;^) are more complete. Indeed, the electron 
positron annihilation R-ratio was calculated long ago |Q in the limits m/ = {f = u,d, s,c, b) 
and rut —>■ oo. The corrections due to the large top-bottom mass splitting have been evaluated 

"'^The part of this work has been done in collaboration with D. E. Soper. 



in 1^ and thus the axial channel, applicable to Tz^hadrons, has also been validated. The 
O(a^) effects of the virtual heavy quark in the decays of the Z-boson have been evaluated in 
IPHI, PT|. The same effect has been studied previously in in the limit ^ oo. The present 
knowledge of the high order QCD and electroweak corrections to Tz^hadrons are summarized 
in the review article [0], providing all essential details. 

The subject of the present work is the correction due to the nonvanishing "light" quark 
masses. The following discussion will be for the quark of flavor / (/ 7^ t) in general. However, 
in fact, only / = 6 is a relevant case and the masses of u, d, s and c quarks can safely be ignored 
at the Z mass scale. Note that the corrections ~ a^m^/Mf for the vector and axial parts of 
the Z bb decay rate were obtained in |T^ and |T^ correspondingly. Those evaluations are 
based on an indirect approach, using the renormalization group constraints and the results of 
earlier calculations of the correlation functions in the vector |T^, |TB[ and scalar channels |]T7 . 
One of the aims of the present work is to obtain the quark mass corrections to T z^qjq^. in 
both channels by a direct calculation and to check the method and the results used in [O, in 
at once. Moreover, it should be stressed that there is a disagreement between the results of 
111 and pf. I 

The quantity Rz is defined as the ratio of the hadronic and electronic Z widths 



^f=u,d,s,c,b ^Z^qjqj 



The partial decay width can be evaluated as the imaginary part 
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where the function 11 is defined through a correlation function of two flavor diagonal quark 
currents 

1 1 dSe''^^ < Tjlix)£iO) >o= 9,Mmf, Q') - Q^Q.n'(m^, Q^). (3) 

Here, is a large (~ — M|) Euclidean momentum. According to the Standard Model, the 
neutral weak current of quark coupled to Z boson is 



/ 



2V2 



i9f qfii.qf + 9f qf7t^l5qf) 



(4) 



where the electroweak vector and axial couplings are defined in a standard way: 



21 f - 4e/ sin^ 9^^, 
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(3) 



^Numerically the disagreement is not large. I n [flSl the result of was confirmed. However, as it will be 
shown in the present paper, the initial result of |_5| turns out to be correct. Also note that unfortunately, in 
lllSi, the result of KM was used. 
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The Il-function may be decomposed into vector and axial parts 

n(m/, = n^(m/, Q^) + n^(m/, Q^). (5) 
The Feynman diagrams that contribute in 11 to O(a^) are shown in Fig. 1. 



Fig.l Feynman diagrams contributing in the correlation function to O(a^). The cut 
diagrams contribute in Tz-,qfgj, at the same order. 

The effects of the last two diagrams in Fig. 1 (so called triangle anomaly type diagrams) 
have been studied in [|, ^ and will not be considered here. In the first diagram, the shaded 
bulb includes any interactions of quarks and gluons (or ghosts) allowed in QCD and the dots 
cover any number of gluon propagators that gives one-, two- and three-loop topologies. The 
crosses denote current vertices corresponding to the vector or the axial parts of the current 

(I)- 

Because the problem scale (~ Mz) is much larger than the quark masses, the following 
expansion is legitimate: 

n^/^(m„m.,g^) = nr(g^) + gin:^(g^)+ ^ 9in:f(g^) + ... (6) 

^ v=u,d,s,c,b ^ 

The last term in the above expansion is due to the certain topological types of three-loop 
diagrams containing virtual fermionic loop. The effects of the virtual top quark in the decays 



of the Z-boson have been studied in |]T0|, |TT|, O and will not be discussed here. The period 



in eq. (^ covers the terms ~ / g^ and higher orders. Those terms at —Q^ = M| for the Z 
decay are heavily suppressed and can safely be ignored. 



The expansion coefficients Ilj can be calculated in a way similar to the one used in for 
the evaluation of the fermionic decay rates of the Higgs boson. In fact, the whole calculational 
procedure can be combined in one equation (in the limit rrit oo) 
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n+k<l 
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s + iO) 



--Ml 



(7) 



where B labels the unrenormalized quantities. Zm and are the MS renormalization con- 
stants of the quark mass and the strong coupling correspondingly and can be found, for 
instance, in |T^. The summation over the virtual quark flavors v = u,d, s, c, b is assumed in 
n*. Note that the introduction of the so called D-function (see, e.g., Q)? which is the 
Il-function differentiated with respect of Q^, is not necessary in this calculation. Note also 
that eq. (|^) does not include important effects from the last two diagrams in Fig. 1 (evaluated 
and the virtual top quark effects (evaluated in [|1^, |T^). 
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In the MS [^] analytical calculations of the one-, two-, and three-loop dimensionally reg- 
ularized Feynman diagrams, the FORM ||2^ program HEPLoops is used. 

For the massless limit coefficients U^^"^ in the expansion (0), the known results are obtained. 



The perturbative expansion of the 



part in the r.h.s. of eq.(g) has the form 



^n:;^(„,(,),Q^)+ i: 

^ v=uA,s,c,b 



i=0 j=0 



7r 



^^^^ ■^{j^'fW'fj^ + ^ij J2 ^lif^))- (8) 



The coefficients Cij are the same in both channels for obvious reasons. Moreover, they get 
nonzero values starting at the three-loop level {i >2). The summation index j runs from zero 
to 1 since, in general, the maximum power of the pole that can be produced by a multiloop 
Feynman diagram is equal to the number of loops. 

The direct computation of all relevant one-, two- and three- loop Feynman diagrams for the 
standard QCD with the SUc(3) gauge group gives in the vector channel: 



d^ 
"oo 



-6; 



d^ 

"lO 



-16, 



-12; 



''20 



19691 124^,^, 1045^,^, 95 ,y 
-^-^C(3) + ^C(5) + iV,-, dl 



253 13 ,y 



Coo — Cij 



0, 



620 



32 



-8C(3), 621 = 622 = 0. 



22 



(9) 
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Note that dYi_^i = Cj^j+i = 0, because the highest poles cancel at each order after the sum- 
mation of Feynman graphs within each gauge invariant set. This is the consequence of the 
conservation of current. The above results fully confirm the findings of |T5[ (see also [^). On 



the other hand, the ({3) coefficient in ^20 disagrees with the incorrect one presented in [|T6|]. 



Fortunately, the numerical difference is small. The error in |T6[ is due to the misprint in the 
C(3) term in the result for the three-loop nonplanar type diagram. 

It can be shown that the vector part of the l.h.s. of eq.(§) is invariant under the renormal- 
ization group transformations and obeys the renormalization group equation 

d d d \ 

'^^'dfJ? ^("^)"«^ ~ lm{as) Yl '^^'q^) ['^^'^tor part of the l.h.s of eq.(8)] = 0, (10) 

where I3{as) and 7m (a^) are the QCD /3-function and the quark mass anomalous dimension 
correspondingly. From eqs.(^|TO|), it is straightforward (similar to ||19|) to get the 0(q;^) 
logarithmic coefficients in the vector channel 

4i = 2(/3o + 7oX + (A + 271)^^1; + 272<o, 

42 = (/3o + 70) [27i4o + (/3o + 27o)rfro] + (A + 27i)7o4o, (H) 
2 

4 = -flM + 7o)(/5o + 27o)4o' c?34 = 0, 
631 = 2{(3q + 7o)e2o, 632 = 633 = 634 = 0. 

The known perturbative coefficients /5„ and 7„ of the QCD /^-function and the quark mass 
anomalous dimension 7^ with the proper normalization factors may be found, e.g., in ||19[. The 
only missing coefficients at the 0(0;^) in eq.(||) for the vector channel are the nonlogarithimc 
terms d^Q and 630. However, these terms have a zero imaginary part and do not contribute in 
the decay rate to O(a^). 

In the axial channel, the direct computation of the relevant one-, two- and three-loop 
Feynman graphs yields: 

6 

doo = - + 6, < = 6; 
< = -^ + - + ^-24C(3), rff\ = 22, 4 = 6; 

. 19 99 /455 , Al 3241 . s 3 , , 

< = 2? - 4? - C(3)) J + — - 387C(3) - 5C(4) + 165C(5) 



/ 1 5 2 857 32^,^, , 



■^in contrary to the previous behef (see, e.g., comments in ||l3, 18 1), that the result of ||l5[ has been corrected 
in 111. 
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where e = {A — D)/2 is the deviation of the dimension of space time from its physical value 4 
within the dimensional regularization Note that the nonlogarithmic terms in df^ contain 



poles, which cannot be removed by the renormalization of the quark mass and the coupling 
and have to be subtracted independently. However, one may not worry about those poles, 
since the imaginary parts of nonlogarithmic terms vanish anyway and do not contribute in 
the decay rate. 

The above mass corrections to the three-loop correlation function of the axial vector quark 
currents are the new results of the present paper. 

One may try to use the renormalization group arguments to obtain the O(a^) logarithmic 
terms similarly to the vector channel ( eq.(0) ). However, to do so, the knowledge of the 
O(a^) anomalous dimension is necessary along the calculated o?^ coefficient. In fact, as it was 



discovered in using the axial Ward identity, this anomalous dimension can be connected 
to the correlation function of the quark scalar currents which, however, is also known only to 

From eqs.@-(|^), (0) and (p!2D, one obtains for the decay rate 



The massless limit coefficients calculated up to 0{al) The coefficients hij read: 



6i; = 12, 6i; = 0, 6fo = -22, 6fi = -12; 
>53 13 

^ ^ y 



&ro = ^-^Ar,, 6^, = 57-2iV,, 6,^, = 0, (14) 



8221 / 1 ^1 \ 1 fi ^7 

h^o = -^+57C(2) + 117C(3)+iV^(— -2C(2)-4C(3)), = -155+-iV^, = --+Nf 

The 0{af) coefficients for the vector part read 

^ 2522 - f C(2) + 5Hc,3) - 5f C(5) - N,{^ - 34C(2, . |^C(3) - ^C(5) 

+iV|(f -HC(2)) (15) 



4505 175 ^^ 13^. ^ 855 ^^^^ 1 . ^ 
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The coefficients c,-i in both channels are 



cij = C2j = 0, C30 = -80 + 60C(3) + Nf 



2? - 5c(3) 

9 3^^ ' 



C31 — C32 — C33 — 0. 



(16) 



The evaluation of 0{a^) coefficients 6^ for the axial part requires the corresponding four-loop 
calculations. The ({2) terms in the above coefficients are due to the imaginary part of the 
term ~ log'^(/i^/s), which appears in the 0{a^) coefficients of the correlation function. 

Taking fj, = Mz, Nf = 5 and recalling the known massless limit coefficients 0, one obtains 
numerically 



8V27T I 



2J[^^ -4e.sin'evK 
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asjMz) 
fas{Mz 



1 + 



1 + 12 



2m'j{Mz 

m){Mz) 
Ml 



\ 



Amj{Mz 



TT 

as{Mz) 

TT 



1.4092 + 104.833 



mjjMz] 
Ml 



-12.805 + 547.879 - 6.12623 



Ml 



E 

v=u,d,s,c,b 



mljMz 
Ml 



+ (2/ 



4m}(Mz)\3/2 



Ml 



+ 



+ 



TT 

as{Mz 



1 - 22 



m){Mz] 
Ml 



TT 

as{Mz) 

TT 



1.4092- 85.7136 



m){Mz 



Ml 



-12.767 + (unknown) 



m){Mz) 
Ml 



6.12623 

v=u,d,s,c,b 



mljMz 
Ml 



;i7) 



where for the Born terms, their well known exact expressions are used. (These terms 
have once again been reevaluated here.) It should be stressed that, in order to obtain a 
complete (up to date) Standard Model expression for the decay rate, the following known 
QCD contributions should also be included: (i) The O(a^) corrections due to the large mass 
splitting within the t-b doublet [§, |[; (ii) The O(a^) effects due to the virtual heavy quark 
0, [n], |12l; (iii) The 0{al) corrections coming from the triangle anomaly type graphs in the 
limit 00 . One also needs to include the electroweak corrections. All those corrections 

can be found in 
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The calculated quark mass corrections to 0{al) and O(a^) gave about 10-20% corrections 
to the corresponding massless results and are of marginal importance for the high precision 
analysis at LEP. It is reasonable to expect that the missing correction in the axial part 

will be of the order similar to the corresponding vector part result. However, the calculated 
mass corrections are important in the low energy analysis, e.g., at PEP and PETRA (or 
B-factory), where the vector part of eq.(17) is relevant. 

For the Z— > bb decay mode, the O(a^) mass corrections agree to the ones obtained in |13, 14 



using an indirect approach, based on the renormalization group arguments and the results of 
|16| , |17|. However, at the 0{al), there is a small disagreement. This, in fact, is due to the 



incorrect numerical coefficient for ({3) term in [^, which was used in [|T3]. Q In the previous 
equations, the strong coupling as{Mz) and the quark mass mf{Mz) are understood as the 
MS quantities renormalized at the Z mass. The relation between the MS running quark mass 
and the pole mass is derived from the on shell results of p5[ (see [|^) 



1 

+ log^ 



/157 13, \, /i^ 



2 



m 



where is the pole mass of the quark, the superscript indicates that the corresponding 
quantity is evaluated for the active flavors of quarks, fi is an arbitrary scale. (In the case 
of Z decay, /i = Mz and A^ = 5.) 



AT/ = 16.00650- 1.04137A^+ - ^ ^(—), 5(mj,m//) 



-1.04137 +^Af!^ 
3 yrrif 



(19) 



and the numerical values for the A at the relevant quark mass ratios are given in |T^. Nu- 
merically, in the case of bb decay mode, Kb ~ 12.5 f\ and the sum over mf drops out in 
eq-dH- 

The calculated mass corrections to the correlation functions are relevant for the hadronic 
decay rates of the r lepton. 



Acknowledgments 

It is a pleasure to thank D. E. Soper for the collaboration, K. G. Chetyrkin, N. G. Desh- 
pande, X. G. He, B. A. Kniehl, M. A. Samuel and D. M. Strom for comments at various stages 
of the work. This work was supported by the U.S. Department of Energy under grant No. 
DE-FG06-85ER-40224. 

"^There is also a misprint in eq.(23) for the general expression for the 0{a^) term in jT^ : the division factor 
92 should be replaced by 96. 

^slightly higher than the one given in p^ . 
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